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Abstract 

We compute the lepton flavour violating couplings of Higgs bosons in the Min- 
imal Supersymmetric Standard Model, and show that they can induce the decays 
fiT at non-negligible rates, for large tan (5 and sizeable smuon-stau mix- 
ing. We also discuss the prospects for detecting such decays at LHC and other colliders, 
as well as the correlation with other flavour violating processes, such as r ^ /X7 and 
r — > 
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1 Introduction 



The recent important indications of neutrino oscillations |JL| reveal that flavour violation also 
occurs in the lepton sector and further motivate the search for alternative signals of lep- 
ton flavour violation (LFV). The Minimal Supersymmetric extension of the Standard Model 
(MSSM) is a natural framework where several such signals could be signiflcant, provided the 
mass matrices of the leptons and of the sleptons are not aligned. Well known examples are 
the LFV radiative decays of charged leptons, fj, — *• 67, r /i7, r — > ej. In this Letter we 
would like to explore another class of such processes, namely the LFV decays of the neutral 
Higgs bosons {h^,H^,A^). An important feature of these decays is that the corresponding 
amplitudes do not vanish in the limit of very heavy superpartners, since the leading con- 
tributions are induced by dimension-four effective operators, at variance with the case of 
radiative decays. 

Related investigations on flavour violating Higgs couplings in the MSSM framework have 
mainly focused on processes with virtual Higgs exchange (see e.g. |21 El IH Ej) and regard 
either quark or lepton flavour violation. The decays of physical Higgs bosons into fermion 
pairs have been explored in the case of quark flavour violation in the MSSM pi 171, whereas in 
the case of lepton flavour violation existing studies |8j |9j have mainly used phenomenological 
parametrizations of the LFV couplings^. 

This Letter is organized as follows. In Section 2 we present the effective LFV Higgs 
couplings in the MSSM framework, focusing on the second and third lepton generations. We 
explicitly compute the one-loop contributions to those couplings and the branching ratios of 
the decays {h^, H^, A^) fir. New results on flavour conserving Higgs couplings are also 
presented. In Section 3 we give a numerical discussion on the LFV Higgs couplings and 
branching ratios, and also discuss the prospects at future colliders. Finally in Section 4 we 
comment on the correlation of the LFV Higgs decays with other LFV processes, such as 
T ^ fi'j and r 3/i, and summarize our results. 

2 Higgs- muon-tau effective interactions 

The MSSM contains two Higgs doublets Hi and H2, with opposite hypercharges. Down- 
type fermions, which only couple to Hi at the tree level, also couple to H2 after the inclusion 
of radiative corrections ^Tj. In particular, for the charged leptons of second and third 
generations the tree-level couplings read as 

C = -Y^Hl^f^i - YrHyr + h.c, (1) 

where Hi is the neutral component of Hi and Yn,Yr are the Yukawa coupling constants^. 
Also the leading effective interactions with H2, which arise once superpartners are integrated 

^An attempt to study LFV Higgs decays in the MSSM can be found in JU]. However, we believe that in 
this work the Higgs couphngs to the sleptons have not been properly identified. 

^We adopt two-component spinor notation, so /i and r {p,'^ and f^) are the left-handed (right-handed) 
components of the muon and tau fields, respectively. Throughout our discussion we assume CP conservation 
and therefore all the dimensionless as well as dimensionful parameters are taken to be real. 
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out, are described by dimension-four operators. These can be either flavour conserving (FC): 



or flavour violating (FV): 

ACfV = -YrALHl^T^^l - YrARH^*lj'T + h.c, 



(2) 



(3) 



where A^, A^, A,- and Al, A/^ are dimensionless functions of the MSSM mass parameters, 
to be described below. In eqs. (j21) and Q we have only retained the dominant terms, 
proportional to Yr, besides the first term in ACpc proportional to F^. In the following we 
are mostly interested in the effects induced by the terms in Q. In the mass-eigenstate basis 
for both leptons and Higgs bosons, the FV couplings read as: 



FV 



y2cos/3 



(Ai r> + Ar /iV) cos(/5 -a)-H^ sin(/5 - a) - iA°] + h.c, (4) 



where tan/3 = {H2) / {H^), a is the mixing angle in the CP-even Higgs sector [V2Re{H^ - 
{H^)) = cos a-h^ sin a, V2Re{H^ - {H^)) = sin a + h° cos a] and A° is the physical 
CP-odd Higgs field. The expression in holds up to 0{Ar tan j3) corrections, which arise 
from eq. and can be O{10%) for large tan/?. For our purposes it is not compelling to 
include and resum such higher-order (tan /3-enhanced) terms. 

The effective couplings (jl)) contribute to LFV low-energy processes, such as the decay 
r — s> 3/i and other ones, through Higgs boson exchange jSl IH E] . We will comment later on 
r — s> 3^. Here we are interested in a more direct implication of those LFV couplings, i.e. the 
decays <I>° fi^r"^ where = h^, It is straightforward to compute the branching 



ratios BR{^^ 



the flavour conserving decays $ 
where the C$ coefficients are: 



and it is convenient to relate them to those of 



r r 



Ch = 


cos(/3 — a) 


2 

, Ch = 


sm{P — a) 




sin a 




cos a 



Ca = 1. 



(5) 



(6) 



Since non-negligible effects can only arise in the large tan/3 limit, in eq. (0) we have approx- 
imated 1/ cos^/3 ~ tan^/9. 

We now present explicit expressions for the quantities Al and A^, i.e. the coefficients 
of the dimension- four operators in (jH)). The relevant one- loop diagrams, which involve the 
exchange of sleptons, gauginos and Higgsinos, are shown in Fig. 1. The diagrammatic compu- 
tation is consistently performed in the gauge symmetry limit, at zero external momentum^. 
In the superfield basis in which the charged lepton mass matrix is diagonal, the mass matrices 

■^In particular, the only Higgs insertion we consider is that explicitly depicted in the diagrams. Fur- 
ther Higgs insertions or momentum dependent effects correspond to higher dimension operators and give 
subleading corrections to $° — > /xr, in the limit of heavy superpartners and large tan (3. 
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Figure 1: Diagrams that contribute to Al [(a),(b), (c),(d)] and to A/j [(e), (f)]. 



of the left-handed and right-handed sleptons read: 
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We are interested in scenarios with large LFV, either in M\ [(LFV)l] or in M\ [(LFV)/?]. 
Large (LFV)l means that fh\^^ is comparable to m\^^ and rh\^^. Similarly, large (LFV)r 
means that ml^^ is comparable to and rh\^^. The flavour states = (z/^, /Ul)"^, L,- = 

{^ti'^l)'^ are related to the M\ eigenstates L2 = (i>2, e^j)"^, L3 = {1)3,6]^^)^ by the relations 
= CLL2 — SLL3, Lt = slL2 + clL^. Analogous relations hold for the right-handed sleptons: 
P'R = crCr^ —srCji^, fji = srCr^ + crCr^, where cr^ and cr^ are the eigenstates of M^. The 
mixing parameters satisfy the following relations: 
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R3 



where rnj and m% 



a = 2, 3) are the eigenvalues of M\ and M% respectively. The other 
relevant parameters are the Bino {B) mass Mi, the Wino (H^", H^^) mass M2 and the fi 
parameter. The latter appears in the Higgsino mass terms -^{H'^Hl - H^H:^) + h.c. and 
in the cubic interaction —YrfiH2*f^fL + h.c. The explicit evaluation of the diagrams gives 
for A^: 
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+4 {l{Ml m?,J - I(M^ m?j3; 
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The function I, which has mass dimension —2, is the standard three-point one-loop integrah 



iix^y^z) = . 

{x-y){z-y){z-x) 



(13) 



Our results for the LFV diagrams in Fig. ^ can be compared with similar ones presented 
in [21 El- However, one notices some differences in those works: i) there LFV effects were 
treated at linear order, through the mass insertion approximation; ii) only LFV in the left- 
handed sleptons was considered, since LFV was related to the seesaw generation of neutrino 
masses; in) the relative signs between the B diagram and gaugino-Higgsino diagrams differ 
from ours. This sign is crucial to correctly determine the interference effects, as we will 
see below. Notice that such a sign discrepancy does not depend on the fact that we use a 
different sign convention for the fi parameter. 

For the sake of completeness we also present the expressions of the FC parameters 
A^, A'^, Ar, which are relevant for establishing the relations between the lepton masses (m^, 
mr) and the corresponding Yukawa couplings. Such quantities are induced by diagrams 
analogous to those in Fig. 1 but with the same flavour in the external fermion lines (either 
muon or tau flavour): 



A, 
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+441(^1, ^L'^L) + 441(^1, ^L'^L) + ^4i(^i,/^',^L) 
+^4i(^i , ^L) - 41(^1 , /^', - 4 W> 
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(16) 



These formulas are quite general as they include possible LFV in the slepton mass matrices. 
By setting cl = cr = 1, sl = sr = one can easily recover for A^ and A,- the corresponding 
cases^ without LFV, whereas A^ vanishes as this term requires both (LFV)l and (LFV)ij. 
Incidentally, notice that A^ in eq. © is multiplied by Yr. Thereby, if (LFV)l and (LFV)r 
are both large, the relation between the muon mass and Yukawa coupling could receive 
large (tan/3 and Yr/Y/j^ enhanced) corrections^, and the ratios BR{^^ fi~^fi~)/BR{^^ — >■ 
r^r^) could differ significantly from the tree level expectation (m^/m^)^. However, having 
simultaneously large (LFV)l and (LFV)^? does not seem very natural if the smallness of 
fn^/mr is related to an underlying supersymmetric flavour symmetry. 



3 Numerical results and implications at colliders 

Now we give some numerical examples to appreciate the size of the effects we are discussing. 
For definiteness, we discuss separately the case of large (LFV)l, with negligible (LFV)r, 
and the complementary case of large (LFV)/?, with negligible (LFV)^. Let us redefine in 
((Tj) m\^^ = m\ and m1^^ = ih\. As a representative case of large (LFV)^, we choose 
m\^^ = m\ and m\^^ = 0.8 -mj^, while ~ 0. Analogously, for the case of large (LFV)r 

we choose m1^^ = m\ and mfj^^ = 0.8 ■ m^, while m\^^ ~ 0. We show the quantity iSOA/^p 
as a function of \^\/mL in Fig. |21 and |50Aijp as a function of \^\/mR in Fig. El for fixed 
values of other mass ratios. We have inserted a factor 50 to make it easier the numerical 
estimate of eq. (0) for the reference case of tan/5 = 50. The curves depicted exhibit a 

^In this limit of vanishing LFV, different expressions for A,- can be found in the hterature | 12II13[ I^. and 
some discrepancies exist among them. Our result for A,- is consistent with that in taking into account 
that we use an opposite sign convention for the parameter and include left-right slepton mixing at linear 
order. To our knowledge no explicit expression for A^ or A^ appears in the literature. In principle A^ can 
be distinct from A,-. 

^In this limit of large (LFV)l and (LFV)^;, analogous enhancement effects also appear in the muon 
magnetic and electric dipole operators, see e.g. For similar enhancement effects in the relation between 
quark masses and Yukawa couplings, see e.g. [7]. 
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Figure 2: The quantity |50A/,p as a function of \^\/rhi^ for mj^^^ = 0.8-m| and four choices of the 
other relevant mass ratios: 1) Mi = M2 = uir = rfiL (sohd line); 2) Mi = mL/2>, M2 = friR = rfiL 
(dotted line); 3) Mi = M2 = rhi, fiiR = mL/2> (dashed line); 4) Mi = M2 = rfiL, rfiR = 3mL (thin 
solid line). 




Figure 3: The quantity iSOA^^p as a function of \fi\/mji, for rh^^^ = 0.8 • and four choices of 
the other relevant mass ratios: 1) Mi = fhi = (solid line); 2) Mi = ttir/S, friL = friR (dotted 
line); 3) Mi = rriR, fiiL = itir/S (dashed line); 4) Mi = rfiR, ifiL = ^rhji (thin solid line). 
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common behaviour^ with respect to the ratio \^\/mL or |/i|/m/j: for each curve there is a 
deep minimum which separates the right-side region, where the pure diagram dominates 
as that mass ratio increases (diagram (a) for (LFV)^ and diagram (e) for (LFV)/? in Fig. 
from the left-side one in which the Higgsino-gaugino diagrams dominate. The deep wells 
for either jA^p or |A/jp are due to the destructive interference of the above mentioned 
diagrams. Notice that the interference would be constructive if the sign of Mi were opposite 
to that of Ma. 

In the case of (LFV)l we can see that values of |50Al p larger than ~ 5 x 10"'^ are achieved 
both in the left and right ranges in Fig. |21 The example with niR = rfiL/S (dashed line) 
provides larger values in the range \fi\/mL > 3 since the pure diagram is further enhanced 
by the smaller rhu- In the case of (LFV)/?, values of |50A^p larger than ~ 5 x 10"'^ can be 
obtained for large values of |/i|/m^ (see Fig. E)). An enhancement appears for = rhji/3 
(dashed line), in analogy to the (LFV)^, example mentioned above. On the other hand, in 
the left-side region the values of |50A/{p are smaller with respect to the analogous ones of 
iSOA^^p. Indeed, in this range |A^p is dominated by the H-B diagram (proportional to g''^), 
while is dominated by the H-W diagrams (proportional to g'^). 

We now make contact with the physical observable, i.e. the i?i?($° — > ^'^t^) in (jSJ, 
and discuss the phenomenological implications. We recall that the Higgs boson masses and 
the angle a in the coefficients C$ are also affected, through radiative corrections, by a set of 
MSSM parameters not involved in the determination of A^, A^j, such as the mass parameters 
of the squark-gluino sector (see e.g. {T5j and references therein). The latter parameters 
indirectly affect also the BR{^^ —>■ r+r~) through radiative corrections to BR{^^ bb) 
(see e.g. |16l I13j ). We do not make a definite choice of those parameters and only outline 
some general features of i?i?($° — > fi~^T~) at large tan/3 and the prospects for these decay 
channels at the Large Hadron Collider (LHC) and other colliders''. It is convenient to 
schematically separate the three Higgs bosons into two groups. The CP-odd and one of the 
CP-even Higgs bosons have about the same mass, non-standard (enhanced) couplings with 
down-type fermions and suppressed couplings with up-type fermions and electroweak gauge 
bosons. These bosons, which are mainly contained in H^, correspond to H'^, A° (/i°, ^4°) 
for m^i > (m^i < m*), where ~ 110 — 130 GeV. The other CP-even Higgs has a 
mass ~ and Standard Model-like couplings with up-type fermions and electroweak gauge 
bosons. It is mainly contained in H2 and corresponds to {H^) for ^ fn* (^A ^ ''^*)- 
Let us discuss the different Higgs bosons, assuming for definiteness tan /5 ~ 50, |50Ap ~ 10~^ 
(A = Ai or A/j) and an integrated luminosity of 100 fb~^ at LHC. 

1. If $° denotes one of the 'non-standard' Higgs bosons, we have (7$ ~ 1 and i?i?($° 
T+T"-) ~ 10~^, so i?i?($° —>■ fi~^T~) ~ 10~^. The main production mechanisms at LHC 
are bottom-loop mediated gluon fusion and associated production with bb, which yield 
cross sections a ~ (10^ 10^,20) pb for ~ (100,200,300) GeV, respectively. The 
corresponding numbers of —>■ fi^T~ events are about (10^, 10'^, 2 ■ 10^). These esti- 
mates do not change much if the bottom Yukawa coupling Yf, is enhanced (suppressed) 

^This behaviour would be visible for all the curves if we had not cut the axes. 

^For recent discussions and references on supersymmetric Higgs phenomenology see also |17| . An uncon- 
ventional scenario has been recently discussed in [TH|. 
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by radiative corrections, since in this case the enhancement (suppression) of a would 
be roughly compensated by the suppression (enhancement) of BR{^^ yU+r^). 

2. If $° denotes the other (more 'Standard Model-like') Higgs boson, the factor C$ • 
i?i?($° T^T^) strongly depends on m^, while the production cross section at LHC, 
which is dominated by top-loop mediated gluon fusion, is a ~ 30 pb. For uia ~ 
100 GeV we may have ■ BR{^'^ r+r") ~ lO'^ and 5i?($o /x+r") ~ 10-^ 
which would imply ~ 300 ^'^t~ events. The number of events is generically smaller 
for large niA since Cq> scales as l/m\, consistently with the expected decoupling of 
LFV effects for such a Higgs boson. However, an enhancement can occur under certain 
conditions. In particular, for a range of itia values the (radiatively corrected) off- 
diagonal element of the Higgs boson mass matrix could be over-suppressed. In this 
case the ^^bb'^, ^^tt"^ couplings would also be suppressed and as a result the number 
of fi~^T~ events could be even C(IO^). 

The above discussion suggests that LHC may offer good chances to detect the decays 
•l'^ — s> fiT, especially in the case of non-standard Higgs bosons. This indication should be 
supported by a detailed study of the background (which is beyond the scope of this work), 
for instance by generalizing the analyses in p. At Tevatron the sensitivity is lower than 
at LHC because both the expected luminosity and the Higgs production cross sections are 
smaller. The number of events would be smaller by a factor 10^ — 10^. A few events 
may be expected also at e~^e~ or fJ,'^fJ,~ future colliders, assuming integrated luminosities of 
500 fb~^ and 1 fb~^, respectively. At a collider an enhancement may occur for the 

non-standard Higgs bosons if radiative corrections strongly suppress Yf,, since in this case 
both the resonant production cross section [a ~ (47r/m^)i?-R($° /i"*"/!")] and the LFV 
branching ratios BR{^^ — >■ fi'^T~) would be enhanced. As a result, for light rriA, hundreds 
of fi~^T~ events could occur. 

Finally we notice that, although we have focused on LFV decays of neutral Higgs bosons, 
also charged Higgs bosons have LFV decays, i.e t^^h and —>■ fi^Ur (and related 

charge conjugated channels). Also these decays are controlled by the parameters and 
Ar, at lowest order in SU{2)\y breaking effects. The FV couplings with the charged Higgs 
bosons emerge by taking into account the SU{2)w completion of eqs. (H)) and 0. It is 
straightforward to find BR{H+ r+u^) = tan^/^ jA^p BR{H+ t+u^) and BR{H+ 
^'^Vr) = tan^ (3 I A/jp BR{H^ t^Vt). However, it is more natural to compare ^Ji^Vr 
with — s> /U+z/^ so that: 

BR{H+ ^ li+Ur) ^ f— ) tan2/5 \^r? BR{H+ /i+z/^). (17) 

For tan^/3 [A^jp ~ 10~^ this would lead to a 30% enhancement in the channel ^ fi'^ + 
missing energy. 
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4 Final remarks and conclusions 



A few comments are in order about possible correlations between the decays $° fir and 
other LFV processes. We have seen that non-negligible rates for $° fir can only be 
obtained for large tan/3 and large LFV. In this limit also the decay rate for r fi'-f, which 
is dominated by diagrams analogous to those of Fig. ^ with an extra photon attached PU] , 
is enhanced and could exceed the experimental limit. However, we recall that the rate of 
r — > /i7 decreases as the superparticle masses increase, whereas the rate of $° — > fir does 
not, since the latter is induced by dimension- four effective operators and only depends on 
mass ratios. Therefore to obtain an adequate suppression of r ^ /i7 the superparticle 
spectrum has to be raised towards the TeV region, although some slepton may be lighter. 
For instance, in the case 1) of (LFV)l shown in Fig. |21 (Mi = M2 = niR = ttil), for 
Ifil/TTiL ~ 1 we obtain iSOA^^p ~ 6 x 10~^. In this particular example the present bound 
BR{t fi'-f) < 6 X 10~^ constrains rfiL > 1.4 TeV for tan/3 = 50, which implies 
niin(rfiL2, mi^3) > 0.6 TeV, ma.x{mL^,rhL:i) > 1-9 TeV and Mi, M2,rhR, \fi\ > 1.4 TeV. 

The decays —>■ fir are also correlated to the decay r 3/i. We recall that the latter 
receives tan /3-enhanced contributions of two types: from dipole LFV operators via photon 
exchange and from the scalar LFV operators ^ via Higgs exchange jHl E] . The dipole 
contribution is directly related to the r ^ fi'y decay rate and is consequently bounded, 
i.e. BR{t ~ 2.3 x lO^^ BR{t fi-f) < 1.4 x 10"^. As for the Higgs-mediated 

contribution, we obtain the following estimate: 

.«(.-3.,.^:o-^(^f?)'(l^)'(^5Mid^ 

Therefore, this contribution can exceed the dipole induced one® and be not far from the 
present bound, BR{t 3/i) < 3.8 x 10~^ j2Tj. Notice that the parameter region in which 
this occurs is also the most favorable one for the observation of the — >■ /xr decays, so an 
interesting correlation emerges. 

Throughout our work we have focused on the second and third generations, implicitly 
assuming that large slepton mixing only appears in that sector. In a scenario in which staus 
are mainly mixed with selectrons rather than with smuons, our discussion and numerical 
estimates concerning — > /ir decays can be directly translated to — er decays, with 
obvious substitutions. The case of large smuon-selectron mixing is somewhat different. Al- 
though the strong constraints from /i — 67 can be satisfied by taking sufficiently heavy 
superparticles, the decays — > fie are generically suppressed by the presence of y^. The 
latter decays could be Yr-enhanced if both (LFV)^ and (LFV)/? were present, and staus were 
mixed with both smuons and selectrons. 

^Here our conclusion is in qualitative agreement with that drawn by [H]. On the other hand, the authors of 
conclude that Higgs-mediated contributions to r ^ 3^ are subleading compared to the photonic penguin 
ones. We believe that this different conclusion is partly due to the fact that in "S" the superparticle masses 
are chosen to lie below the TeV scale, so sizeable values for the LFV Higgs couplings are prevented by the 
r — > /X7 constraint. 
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In summary, we have studied the LFV couphngs of Higgs bosons in a general MSSM 
framework, allowing for generic LFV entries in the slepton mass matrices, but without in- 
voking any specific mechanism to generate them. We have computed the branching ratios 
of $° — > fiT decays, which depend on ratios of MSSM mass parameters, and increase for 
increasing tan/3 and LFV. Although cancellations can occur in some regions of parameter 
space, O{10~^) values are achievable, and they are compatible with the bounds on t ^ fi'j 
for a superparticle spectrum in the TeV range. If the Higgs spectrum is relatively light 
{niA ^ 300 GeV), our results indicate that future colliders (in particular LHC) may be able 
to detect the decays $° — *• /xr, especially in the case of non-standard Higgs bosons. More- 
over, the detection of these decays is closely correlated with that of r — > 3fi, which may be 
observed in the near future. 
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